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ABSTRACT 

We have obtained Integral Field Spectroscopy for 23 spiral bulges using INTEGRAL on the William 
Herschel Telescope and SPIRAL on the Anglo-Australian Telescope. This is the first 2D survey 
directed solely at the bulges of spiral galaxies. Eleven galaxies of the sample do not have previous 
measurements of the stellar velocity dispersion (<7„). These data are designed to complement our Space 
Telescope Imaging Spectrograph program for estimating black hole masses in the range 10 6 — 10 8 Mq 
using gas kinematics from nucleated disks. These observations will serve to derive the stellar dynamical 
bulge properties using the traditional Mgb and Call triplets. We use both Cross Correlation and 
Maximum Penalized Likelihood to determine projected c* in these systems and present radial velocity 
fields, major axis rotation curves, curves of growth and c* fields. Using the Cross Correlation to 
extract the low order 2D stellar dynamics we generally see coherent radial rotation and irregular 
velocity dispersion fields suggesting that cr, is a non-trivial parameter to estimate. 
Subject headings: integral field spectroscopy - galaxies: spiral - galaxies: velocity distributions 
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1. INTRODUCTION 

The increasing number of super-massive black hole 
(SMBH) candidates discovered in the mid 1990s led 
to the investigation of SMBH demographics. It was 
found that various properties of the host bulge correlated 
with the SMBH mass, e . g. bu l ge luminosity and mas s 
ijKormendv fc Richstond 119951 IMagorrian et alj 11998ft . 
stellar velocity dispersion - the M. — a* relat ion - 
l)Ferrarese fc MerrittJ pOOOt iGebhardt et all fcOOO) and 
concentration index ( Gra ham et al.l 2003). Scatter in 
these relations are comparable if enough care is take n 
in deriving the bulge properties (Mar coni 8z Hunj l2003'). 
Such correlations have fundamental implications for both 
SMBH and host galaxy formation and evolution. These 
"secondary" methods also present a valuable practical 
tool to estimate SMBH mass. The "primary" methods, 
used to originally define the relations, determine SMBH 
mass by considering the direct gravitational effects of 
SMBHs. However, each primary method has a drawback: 
proper motion studies can only be carried out in the 
Milky Way; H2O masers need to be orientated correctly; 
gas kinematics may be influenced by non-gravitational 
effects and stellar dynamics must be derived from rela- 
tively low surface brightnesses. The advantages of know- 
ing the exact form of the SMBH - host relations then 
become obvious if trying to sample a large number of 
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SMBH masses across a wide dichotomy of host types and 
redshifts. 

There still exists a fundamental bias toward large 
(~ 1O 9 M0) SMBH mass measurements in the relations 
however. This stems from the relative ease of mass 
measurements in giant ellipticals: it is easier to spot a 
bigger black hole. The low mass end of the relations 
(10 6 M© < M. < 10 8 M o ), which essentially represents 
late type spirals, has yet to be fully populated. This is 
unsurprising when considering the difficulties and com- 
plexities involved with observing and modelling such sys- 
tems. This swing toward high mass SMBHs has led some 
to dis cuss whether these relations are indeed linear ijLaorl 

Hoi. 

If we are to effectively extrapolate these secondary re- 
lations to make conclusions about formation and evolu- 
tionary scenarios, and indeed use them as SMBH mass 
estimators, we need to quell any doubts as to their na- 
ture. To this end we have completed a Space Telescope 
Imaging Spectrograph (STIS) survey of 54 local spirals 
(Hug hes et aDl2003t IScarlata et alJl2004|) using the Hub- 
ble Space Telescope (HST). This data has been used to 
make primary SMBH mass estimates by modelling the 
gravitational potential as represented by the central gas 
kinematics and surface brig htness profiles ( M arconi et alJ 
120031 [Atkinson et al J 12005). However, to populate the 
low mass end of the secondary estimators we need to 
establish the bulge properties. Consequently we have 
initiated integral field spectroscopic observations of the 
original STIS sample. Integral Field Spectroscopy (IFS) 
is a reliable and highly flexible tool capable of quickly 
gathering a large amount of dynamical information. It 
allows the observation of the full 2D Line of Sight Veloc- 
ity Distributions (LOSVDs). 

As there are few published 2D kinematical maps of spi- 
ral bulges we present here an atlas of our data so far in- 
cluding radial velocity fields, major axis rotation curves, 
stellar velocity dispersion (cr*) fields and aperture inte- 
grated velocity dispersions (curves of growth). In §0we 
outline the observations and basic data reduction, while 
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Fig. 1. — Demonstrating the Xcor fits to the data, [top row] The left panel shows the cross correlation function derived from a S/N=23 
spectrum in NGC 4041, while the right panel details the function in a S/N=8 case, [bottom row] The bold lines show the fits using the 
derived parameters to the observed spectra. The lower sections present the residuals of the fits. 



§|3]details how we extracted the 2D LOSVDs. The data 
is presented in § 0] and briefly summarized and discussed 
in §0 

2. OBSERVATIONS AND DATA REDUCTION 

Our observations form the first IFS survey specifically 
applied to a sample of spiral bulges. They positively com- 
plement the SAURON IFS project whic h is directed at 
72 ta r gets including a numb er of spirals l(de Zeeuw et al.l 
l2002t lEmsellem et alJ 120041. From the nature of the 
ill, — a. t relation, and the general character of spiral 
galaxies, we required instruments capable of resolving 
<7* well below 100 km s . For our northern hemisphere 
targets this lead us to use the fibre fed Integral Field Unit 
(IFU) INTEGRAL, on the 4.2m WHT, with the 5 th or- 
der Echelle grating (3? « 6000). Whilst th e Call triplet is 
the preferred tracer of stellar dynamics i|Dresslerlll984l) 
this grating could only be centered on the Mgb triplet 
(~ 5200 A). Observing conditions ensured all nights were 
non-photometric and forced the use of the SB3 fibre bun- 
dle. This bundle consists of a main 33'.'6 x 29'.'4 array 
with 115 circular apertures 9 packed according to Kepler's 
Conjecture. The 2'.'7 diameter of these apertures ensured 
that the PSF due to seeing (which was typically 1'.'5) was 
under-sampled. The main array is surrounded by a con- 
centric circle (090") of 20 apertures designed to sample 
the sky. For our southern hemisphere targets we used 
the fibre fed IFU SPIRAL, on the 3.9m A AT, with the 
R1200 grating (JR. « 7600) centered on the Call triplet at 
~ 8660A. Ordinarily SPIRAL consists of 512 apertures 
but effective sky reduction can be achieved by masking 
half of the array and using SPIRAL in "nod and shuffie" 
(NS) mode. NS masking of the MITLL2A detector meant 
that 1 row and 1 column of the array fell off the chip. 
This reduced the number of usable apertures ((7)0"7) to 
225 over a 11" x 9"8 area. 

During our observing runs priority was given to targets 

9 Throughout "apertures" means the individual fibres, lenslets 
or spaxels referred to in other IFS studies. 



for w hich we have successful STIS spectra l|Hughes et alJ 
2003). To date we have completed IFS observations on 
23 of the original 54 spiral bulges. Calibration frames 
including flat-fields, tungsten lamp flats, dark frames, 
bias frames and Cu-Ne/Cu-Ar arc frames, along with 
at least 3 late type (G, K, M) spectroscopic standards, 
were observed throughout each night. At least 3 tar- 
get exposures were made in order to facilitate the easy 
removal of cosmic rays. We made several offset sky ob- 
servations, especially when the overheads associated with 
NS observations lead us to abandon its use in some cases. 
The SPIRAL data was reduced using the adapted auto- 
matic 2df software developed by Jeremy Bailey of the 
Anglo- Australian Observatory. INTEGRAL data reduc- 
tion was carried out with standard IRAF 10 routines, and 
made use of the integral package developed by Robert 
Greimel of the Isaac Newton Group. Detector positions 
of the individual aperture spectra were traced by fitting 
a polynomial across the dispersion axis of a high signal 
lamp flat. This enabled scattered light to be identified 
and removed from the science frames. The trace of the 
spectra was then used as a master template in extracting 
the subsequent spectra. Lamp flats were again used to 
make throughput corrections before carrying out wave- 
length calibrations. Sky subtraction was performed us- 
ing NS or by subtracting a collapsed off-set sky frame 
or summed sky apertures. All objects had the contin- 
uum fitted to each individual spectra. This level was 
then used together with the aperture spacial positions 
to create 2D intensity maps (essentially very low spacial 
resolution images of the targets). 

3. EXTRACTING LOSVDS 

Since the early 1970s there have been a multitude of 
techniques employed to e xtract dynami c al information 
from digitized spectra l|Simkinl 119741 ISargent et alJ 

10 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under co-operative agreement 
with the National Science foundation. 
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Fig. 2. — Identifying discordant data. In all cases the solid line marks the median value and the dashed lines show the upper and lower 
quartiles. The dotted lines mark the 2i ± m boundaries, where m is the median and i is the inter-quartile range. Data falling outside of 
the dotted lines were marked as discordant, [left column] NGC 4212. The top panel presents the radial velocities and the bottom panel 
shows the velocity dispersions in this low discordant example, [right column] Same as the left column but for NGC 3021; a high discordant 
example. 

the Xcor. Quadratic weights are applied to deviations in 
the continuum subtracted spectra to ensure that stronger 
lines are more heavily weighted. To limit the size of sec- 
ondary peaks in the Xcor function a cosine bell filter is 
applied. 

Errors in the position and width of the Xcor are de- 
fined by the "r" statistic which is calculated by rvsao and 
is used as a quality measure of the Xcor. If we assume 
that a perfect Xcor is proportional to the broadening 
function convolved with the Auto Correlation Function 
(ACF, i.e., the Xcor of a spectrum with itself) then we 
expect to see symmetry about the radial velocity. This is 
because both the broadening function and the ACF are 
themselves symmetrical. The constant of proportionality 
between the perfect correlation and the convolved broad- 
ening function will be defined by the height of the central 
peak. The height of the peak, normalized by the anti- 
symmetric components in the Xcor, defines r and is used 
as the aforementioned quality measure of values derived 
from the Xcor itself. 

However, these error estimates are insignificant when 
compared to those introduced by the possibility of tem- 
plate mismatch, low signal spectra and large velocity dis- 
persions. The advantages of observing a template library 
now becomes apparent. For each aperture, after correct- 
ing the Xcor radial velocities for barycentric and stellar 
radial velocities, we used the r weighted average value 
between each template for the velocity and width mea- 
sured in that aperture. The standard deviation around 
these weighted values was then used as a separate es- 
timate in the error of each measurement. More weight 
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In general each method 
makes the fundamental assumption that LOSVDs can 
be recovered by convolving a broadening function to 
a shifted template stellar spectra of the same mean 
spectral type. We have attempted to use each of these 
methods on our sample, but only a select few are suited 
to the nature of our data (essentially ~ 3500 S/Nw 10 
spectra). Typical Fourier and pixel fitting techniques 
fail at the levels of S/N present here and require a 
large amount of time to complete due to the nature 
of filters and number of possible solutions. We have 
seen t hat the technique of C ross Correlation, hereafter 
Xcor (|Tonrv fe Davies 1119791 TD79), performed well in 
extracting accurate 1 st and 2 nd order LOSVD moments 
from our spectra and is particularly suited to large 
datasets. 

We performed Xcor on each spectrum for each ob- 
served template by following TD79 as implemented in 
the iraf. rvsao package. As these reduction steps are con- 
sidered standard and well established we refer readers 
wishing to see the detailed treat ment of spectra throug h 
the cross correlation procedure to lKurtz fc Minkl (1998). 
In brief, the position of the Xcor peak determines the ra- 
dial velocity shift, whilst the width of the peak is related 
to the projected stellar velocity dispersion (<r*). Each 
spectrum has any residual emission features suppressed 
and is continuum subtracted, apodized, zero-padded and 
binned linearly in log space before being passed through 
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TABLE 1 
Main Results for the Sample 
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Note. — "% DC" refers to the percentage of discordant data. "<r cc " indicates extended Xcor values, 
"<t cc " central Xcor, "<r mc " extended MPL and "cr mc " central MPL. "/13" and "/14" are for the higher LOSVD 
moments in the central apertures. 



was placed on higher r values as defined in the following 
equation: 



where <pp is the final velocity or width and 4>x..3 and 
ri..3 represent the values derived using template 1,2 or 
3. Any of the weighted values that fell outside the equiv- 
alent of a 2.7 <j clip (defined by the interquartile range 
of the data from the entire array) were flagged as "dis- 
cordant" , the percentage of discordant events per galaxy 
then being used as a quality measure of the IFS data. 

As with the intensity maps, the full 2D kinematics were 
recovered by mapping the derived velocities to the rela- 
tive spacial positions of each aperture on the IFUs. 

We have carried out extensive testing of the Xcor pro- 
cedure by creating a composite spectra from the observed 
templates, generating simulated radial velocity and ve- 
locity dispersion fields, introducing a surface brightness 
profile and adding Poisson noise. As the width of the 
Xcor peak does not give a direct measure of cx* we tested 
three independent 0% recovery methods: HWHM correc- 
tion (er\/2m2); ACF correction (c 2 = /1 2 — 2r 2 , where /x 
is the width of the peak and r is given by the ACF) and 
a calibration function correction. The calibration func- 
tion, which was determined by fitting a polynomial to 
the known and measured values of cr*, proved to be far 
superior in terms of x 2 . As some of our sample are known 



to contain starbursts and AGN we also tested the effects 
that such objects may have on the Xcor. To the compos- 
ite template spectrum we added increasing proportions 
of early type spectra and a Seyfert 1 spectrum. In both 
cases we found that 90% contamination was needed in 
order to significantly affect the r statistic (i.e., take it 
below 3). Similar results were also found when introduc- 
ing large velocity shifts to spectra (~ 3000 km s" 1 ). 

Examples of the Xcor are presented in Figures ^ and [21 
Figure ^ shows the cross correlation functions from two 
apertures in NGC 4041 and overlays of the derived fits 
onto observed spectra. In both the high signal (S/N=23) 
and low signal (S/N=8) cases well defined cross correla- 
tion functions can be seen giving confident results. These 
Figures also demonstrate the nature of the r statistic (the 
amount of anti-symmetry in the cross correlation func- 
tion). Figure |2] demonstrates how discordant data was 
identified and flagged. The two cases presented show 
both low and high discordance observations. The de- 
rived radial velocity and velocity dispersion data from 
low discordance and high discordance cases (NGC 4212 
and NGC 3021 respectively) were placed into discrete 
bins in velocity space. The distribution of data is shown 
by plotting the number of measurements to fall within a 
specific bin. 

In spectra down to S/Nss 5 we are confident that the 
Xcor procedure will return accurate values for the 1 st 
and 2 nd order LOSVD moments. However, in many cases 
the central apertures of the IFU arrays contain S/N > 15. 
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Fig. 4.— Same as Figure|3]but for NGC 0613. 



Here we can attempt to recover higher order moments (/13 
and /14) which harbor further dynamical information on 
the observed system, i.e., signatures of SMBHs. Perform- 
ing LOSVD extraction via an independent technique also 
provides a valuable verification of previous resul ts. The 
Maxi mum Penalized Likelihood, hereafter MPL l)MerrittJ 
[1993), is a non-parametric method in that it places no 
constraints on the shape of the LOSVD and allows the 
recovery of higher order moments {ho. ..he). It is ideally 
suited to low S/N as the penalty function introduces sig- 
nificantly less bias than Fourier methods when smoothing 
noise characteristics. We carried out similar tests as for 
the Xcor and again found that 90% starburst/AGN con- 
tamination and shifts of > 3000 km s _1 were needed to 
significantly affect the extracted moments. In order to 
determine the most appropriate value for the smoothing 
parameter (a) we re-performed the MPL for increasing 
values of a until cusps in the higher moments were seen. 
A bootstrap, using 100 realizations, was used in order 
to determine the 68% confidence levels in the derived 
Hermite polynomials. 
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Fig. 5.— Same as FigureEJbut for NGC 1255. 
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Fig. 6.— Same as FigureEJbut for NGC 1300. 

4. RESULTS, VELOCITY FIELDS AND COMMENTS FOR 
INDIVIDUAL GALAXIES 

The main results for our sample are presented in Ta- 
ble n We have provided measures of cr* taken from the 
central aperture of each array (determined by the peak 
value of cr* closest to the photometric center of the target 
- "central" apertures) and also cr* measured from spectra 
collapsed across the entire arrays ( "extended" apertures). 
We also present, for each galaxy, intensity maps derived 
from the continuum level of each spectra, the "curve of 
growth" derived by determining a„ through progressively 
larger collapsed apertures, the radial velocity field and 
the velocity dispersion field. In all cases North is up and 
East is to the left. In the cases where there is evidence 
of a net rotation the radial velocities along the major 
axes have also been presented. We continue by making 
comments on each galaxy. 

4.1. NGC 0289 

Figure The data are noisy and almost 1/3 is ex- 
cluded as discordant. As a consequence nothing can be 
said about the 2D dynamics. These are the first measures 
of cr, in NGC 289. 
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Fig. 7.— Same as Figure|3]but for NGC 1832. 
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Fig. 8.— Same as FigureEJbut for NGC 2748. 

4.2. NGC 0613 

Figures Bland II Uf a*). Clearly defined rotation despite a 
high level of discordance. The rj* field is patchy but a sig- 
nificant difference between extended and central values 
of cr* is noted. There are also significant values of /13 and 
hi corresponding to a central asymmetric non-Gaussian 
LOSVD. 

4.3. NGC 1255 

Figure |3J S /N is low in all spectra but the level of 
discordance is also low. No clear rotation can be seen 
and the rr* field is patchy. These are the first measures 
of cr, in NGC 1255. 

4.4. NGC 1300 

Figures 151 and rTUT b). A low level of discordance and 
good S/N has produced a velocity field showing rotation 
and a clear centrally peaked <r* field. Color maps suggest 
possi ble obscuration from spiral structures l)Hughes et al.l 
2003). 

4.5. NGC 1832 

Figure [7| The levels of S/N are low and the radial ve- 
locity and <7* fields are noisy. There is a significant differ- 
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Fig. 9.— Same as FigureEJbut for NGC 2903. 
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Fig. 10. — Radial velocity curves taken along the major axes 
of each galaxy. Numbers in square brackets refer to the position 
angle in degrees (east from north) along which the curves have been 
extracted, (a) NGC 613 [135]. (b) NGC 1300 [90]. (c) NGC 2748 
[30]. (d) NGC 2903 [150]. 

ence between central and extended values of rj„. These 
are the first measures of ex* in NGC 1832. 

4.6. NGC 2748 

Figures IS1 and HUT c). Typical S/N and discordance. 
Clear rotation and patchy <r*. NGC 2748 has an in- 
clination of ~ 70° which could account for LOSVD 
obscuration especia lly in the presence of dust lanes 
(Hughes et al. 2003). These are the first measures of 0% 
in NGC 2748 and are constant with aperture radius. 

4.7. NGC 2903 

Figures |H1 and HOf d). High level of discordance 
but S/N is good. Clear rotation and an irregular 
cr* field. MPL did not converge in the central aper- 
ture. A known star-burst with a circ umnuclear ring 
onso-Herrero. Ryder fc Kn apcii 2001|). 

4.8. NGC 2964 

Figures ITTI and UTT a) . Good S/N with clear rotation 
and a centrally peaked rx* field. Slight central peak 
in the curve of gr owth. Significant amounts of dust 
lIHughes et aLlEffll . 
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Fig. 11 — Same as FigureEJbut for NGC 2964. 



Fig. 13.— Same as FigureEJbut for NGC 3162. 
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Fig. 12.— Same as Figure|3]but for NGC 3021. 

4.9. NGC 3021 

Figures ITU and IT7T b - ) . Significant discordance and low 
S/N. Clear rotation and a peaked ct* field. Significant 
difference between extended and central apertures. 

4.10. NGC 3162 

Figures IT31 and fPTT c). Low S/N and typical discor- 
dance. Clear rotation and irregular <r* field. Signifi- 
cant difference between extended and central apertures. 
There are central asymmetries in the LOSVD. These are 
the first measures of <r* in NGC 3162. 

4.11. NGC 3310 

Figures d and fTTTd). Good S/N and typical dis- 
cordance. Clear rotation, patchy <r* field and a clear 
difference between central and extended <r*. Possible 
high velocity co mponents in central L OSVD. Significant 
amounts of dust l)Hughes et alJl2003|) . 

4.12. NGC 3949 

Figures IT51 and l2*UT a) . High discordance and fair S/N. 
Clear rotation and irregular er* field. Constant curve of 
growth. Significant high velo city components in the cen- 
tral LOSVD. Thin dust lane (|Hughes et alJl2003|) . These 
are the first measures of cr» in NGC 3949. 
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Fig. 14.— Same as FigureEJbut for NGC 3310. 

4.13. NGC 3953 

Figures [TBI and l2"UT b') . Low discordance and good S/N. 
Centrally peaked er* field and clear rotation. Large cen- 
tral curve of growth peak. There are no previous mea- 
sures of cr* for NGC 3953. 

4.14. NGC 4041 

Figures ^| and I2l)f c). Good S/N and low discor- 
dance. Fairly constant curve of growth, regular rota- 
tion and a centrally p eaked ct* field. Nuclear star cluster 
ijMarconi et al.ll2003(l . There are no previous measures 
of ct* for NGC 4041. 

4.15. NGC 4051 

Figures IT^l and l2*UT d') . Significant discordance and low 
S/N. The large peak signal is due to the AGN. Regular 
rotation, patchy cr* field and a steep curve of growth. 
There are large values for both /13 and /14, the signature 
of a large asymmetry and super-Gaussian wings in the 
LOSVD. 

4.16. NGC 4088 

Figures l2"71 and ISST a) . Low S/N and high discordance. 
Clear rotation and a patchy er* field. A significant cen- 
tral peak in the curve of growth. The inclination of 
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Fig. 16.— Same as FigureEJbut for NGC 3953. FIG. 19.— Same as FigureEJbut for NGC 4051. 




NGC 4088 (~ 71°) suggests there may be some obscura- 
tion from the disk. These are the first measures of er* in 
NGC 4088. 



NGC 4212 

Remarkably low discordance 
and good signal. Centrally peaked cr* field, clear regular 



4.17. 

Figures HI and H^b). 



rotation and a significant peak in the curve of growth. 
There are no previous measures of cr* for NGC 4212. 

4.18. NGC 4258 

Figures EH and |2H C ) ■ Typical discordance and excel- 
lent S/N. Clearly defined rotation and an irregular cr* 
field. NGC 4258 is inclined at ~ 72° and exhibits a dust 
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Fig. 23.— Same as Figure|3]but for NGC 4258. 

lane l)Hughes et alJl2003l) . The curve of growth shows a 
steep drop toward the center. 

4.19. NGC 4303 

Figures [31 and I2tjf d~) . Good S /N and typical discor- 
dance. There is clear offset rotation and a patchy cr* field. 
Significant difference between extended and central cr*. 




Distance ~rorr Centre (") Distance from Cer ,_ .re 



Fig. 26.— Same as Figure ESI but for (a) NGC 4088 [30]. (b) 
NGC 4212 [90]. (c) NGC 4258 [120]. (d) NGC 4303 [120]. 

Nuclear spiral structure ijHuehes et al.l l2003'i . 

4.20. NGC 4321 

Figures I^SJand |20{ a). Low discordance and good S/N. 
Clear peak in <r* offset from the photometric center. 
Clear regular rotation and a centrally peaked curve of 
growth. 
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Fig. 27.— Same as Figure|3]but for NGC 4536. FlG - 29 -— Same as Figure|3]but for NGC 5055. 
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Fig. 28.— Same as FigureEJbut for NGC 5005. 

4.21. NGC 4536 

Figures I2"7land l3*07 b). Significant discordance and very 
low S/N. Steep rotation and centrally peaked a* field. 
Steep centrally peaked curve of gro wth. High inclinatio n 
(~ 69°) and obscuring dust present ijHuehes et alfc OOSI. 

4.22. NGC 5005 

Figures |2HI an d !3Uf cL Low discordance and excel- 
lent S/N. Clear, well defined rotation and complex cr* 
field. Highly variable curve of growth. Large incl ination 
(~ 67 °) and complicated dust morphology ( Hugh es' et al.l 
I200H) . These are the first measures of cr* in NGC 5005. 

4.23. NGC 5055 

Figures l2"51 and l3"07 d). Low discordance and excellent 
S/N. Constant curve of growth and clear, well defined 
rotation. Patchy cr* field. 

5. SUMMARY AND DISCUSSION 

In order to populate the low mass regimes of the SMBH 
- host relations we have completed a STIS survey of 
54 bulges in spiral galaxies from which to derive the 
SMBH mass from ga s kinematics lIMarconi et al.ll2003t 
lAtkinson et al.ll2005j) . To place these points on the re- 
lations we need to determine the dynamical properties 
of the host bulges. To this end we have initiated the 
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FlG. 30.— Same as Figure llOl but for (a) NGC 4321 [120]. (b) 
NGC 4536 [120]. (c) NGC 5005 [30]. (d) NGC 5055 [90]. 

first IFS survey specifically directed toward spiral bulges. 
We have presented an atlas of 23 spiral bulges, eleven of 
which do not have previous estimates of cr* . 

The properties of the LOSVDs were derived using the 
methods of Cross Correlation (Xcor) and Maximum Pe- 
nalized Likelihood. Each method was thoroughly tested 
against a simulated system and are robust to spectral 
contamination from AGN and starbursts. 

Where the S/N has allowed we see that each sys- 
tem shows clear coherent radial rotation and patchy 
irregular velocity dispersion fields. This type of be- 
haviour is displayed in a number of systems observed 
by SAURON, e.g. NGC 4382 (S 0*(s)pec), NGC 4473 
(E5) and NGC 4526 (SABO*(s)) (|de Zeeuw et alJl200l 
lEmsellem et al.ll2004|) . We also note that the derived 
curves of growth, i.e. the variation of cr* with aperture 
size, show large variations of cr* in some cases. Generally 
we see offsets between the photometric centers of these 
systems and their peak values of cr*. We are confident 
that this observation is also robust in the cases where 
a large amount of data has been excluded due to low 
S/N or template mismatching. This discordant data is 
mainly confined to the outer areas of the Integral Field 
Spectroscopy arrays where the S/N is at its lowest. As 
a consequence data at the centers of the arrays are still 
of high quality (S/N > 14.5). There is also a correlation 
between the complexity of cr* fields and the inclination of 
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the system, especially where it is known that dust obscu- 
ration is significant, i.e., in high inclination cases (> 70°; 
NGC 2748, NGC 4088, NGC 4258, NGC5005) the larger 
values of a* are confined to a band coincident with the 
photometric major axis. These observations suggest that 
a* - the projected stellar velocity dispersion - is a non- 
trivial parameter to derive consistently from system to 
system. 

With these considerations in mind we are left to con- 
template the errors in SMBH mass estimates derived 
from the M. — a* relation. The data presented here sug- 
gests that values of cr* taken from the literature, where 
little heed has been taken of aperture sizes and the effects 
of rotation, may demonstrate variations due to their het- 
erogeneous nature. An attempt to quantify these errors, 
with respect to aperture size, can be made by calculat- 
ing 5a*, i.e. the mean difference between the extended 
(cr co ) and central (cr cc ) values of <t* as determined from 
Xcor. Indeed, if we are to assume that a co and a cc are 
a good representation of er„ within R e and R e /8 respec- 
tively (R e being the effective or half-light radius), then 
5a* may also go some way toward explaining the scatter 
in the M, — a* relation, especially in the spiral region 
where the most weighting is given to the slope (a). 

From the data presented in Table Q we find 5a* 
20 km s _1 . This figure implies that the slope of the 
Mm — a* relation will be affected by the aperture size 
used to measure a*. It then becomes clear that a 



formal unambiguous definition of a* must be made, 
undoubtedly through the use of Integral Field Spec- 
troscopy where the exact behaviour of a* can be fol- 
lowed across two dimensions, before a confident esti- 
mate of the slope and scatter in the M. — a* relation 
can be made. It must also be noted that carrying this 
value of 5a* through th e two forms of the M, — a * re- 
lation, i. e., a = 4.72 llMerritt fc Ferraresd l200l|) and 
a = 4.02 ijTremaine et alJl2002|) using 100 - 120 km s -1 , 
leads to SMBH mass ranges of 0.49 - 1.17 x 10 7 M Q and 
0.83 — 1.73 x 10 7 M Q respectively. Therefore, in future, 
authors wishing to use the M, — a* relation to estimate 
the masses of SMBHs from "simple" a* measurements 
must pay close attention to the processes used in order 
to make that estimate of a* in the first place. 
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